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Abstract 

We introduce a new notion of viscosity solutions for a class of very singular nonlinear 
parabolic problems of non-divergence form in a periodic domain of arbitrary dimension, 
whose diffusion on flat parts with zero slope is so strong that it becomes a nonlocal quantity. 
The problems include the classical total variation flow and a motion of a surface by a 
crystalline mean curvature. We establish a comparison principle, the stability under approxi- 
mation by regularized parabolic problems, and an existence theorem for general continuous 
initial data. 



1 Introduction 

We are interested in the well-posedness of the initial value problem for a class of very singular 
parabolic equations. They are characterized by the strong diffusive effect on flat parts of the 
solution with zero slope. In fact, the diffusion is so strong that the equation becomes nonlocal 
and it is no longer a classical partial differential equation. The problems with such structure have 
attracted considerable attention due to their importance in image processing and crystal growth 
modeling. A typical example of a problem of divergence type is the standard total variation flow 
equation. Its well-posedness is not trivial but the theory of monotone operators yields the unique 
solvability of the initial value problem [[TJ. However, this theory does not apply if the problem 
lacks a divergence structure and the unique solvability of the initial value problem has been a 
long-standing open problem except in the case of one space variable for which the theory of 
viscosity solutions is extendable [24|26|. Because of the ubiquity of such higher dimensional, 



non-divergence form problems in crystal growth modeling, it is desirable to find a suitable 
generalized notion of solutions and establish the unique solvability of the initial value problem. 

Our goal in this paper is to solve this long-standing open problem in a periodic domain of 
arbitrary dimension n > 1 by introducing a new notion of viscosity solutions. In full generality, 
we consider a generalized total variation flow equation for a function u(x, t) on the n-dimensional 
torus Q. = V := R"/Z": 

Mf + FlVMiv&)) = in£:=T"x(0,r), (1.1) 
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with the initial condition 



u\ t=0 = u onT", (1.2) 

where u, = § , Vw = (j^, . . . , and divz = g- + • ■ • + ||. We assume that F : R" x R -» R 
is a continuous function, non-increasing in the second variable, i.e., the operator is degenerate 
elliptic: 

F(p,Z)<F(p,T]) for^,^eR,^>77,peR". (1.3) 

The initial data w is assumed to be a continuous (periodic) function uq e C(T). We identify 
functions from C(T") with Z" -periodic functions in C(R"). 

The prototypical example of ( |1.1| ), with many applications in image processing, is the 
standard total variation flow 

/ Vw \ 

u, = div — , (1.4) 



JVh, 

which is formally the L 2 -gradient flow of the total variation energy 

E(iff) = I \Vif/\ dx. 
Jn 

The gradient flow structure can be exploited to employ the abstract theory of monotone operators 
and obtain well-posedness for general data in Hilbert space L 2 ; we refer the reader to nTJ for 
a detailed exposition and references. This interpretation will serve as the basis of evaluating 
div(Vw/ |V«|). 

The motion of an interfacial surface by a crystalline mean curvature is an example of a 
problem in non-divergence form, 

u t = 7l + |Vw| 2 fdiv(^J + cJ, (1.5) 

where c is a constant forcing. This is the graph formulation of the anisotropic mean curvature 
flow of a surface T t = {(x, u(x, t)) : x e R"} c R" x R of codimension one, given as the graph of 
function u, moving with velocity 

Vy = K y + C (L6) 

in the direction of the unit normal vector v = (-Vw, 1)/ yjl + |V«| 2 . The quantity k 7 represents 
the crystalline mean curvature which is the first variation of the surface energy 



I 



y(v) dfW, 



see [TT}36|, where "H" denotes the rc-dimensional Hausdorff measure. The Finsler metric y(x, u) 
has the particular form 

y{x, u) = \x\ + \u\, x e R", u e R, 
and the Wulff shape, which represents the surface of constant curvature K y , is the unit cylinder 

W y = {(*,«) : max(|jc| , \u\) < 1} c R" +1 . 
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We mention that y is not purely crystalline in the traditional sense for n > 2, because usually 
only Finsler metrics with a polytope as their Wulff shape are considered crystalline; see 
for an introduction to the topic. We observe that on the flat parts of u with zero slope, k 7 and 
div(Vw/ |Vw|) coincide. Motivated by this example, we shall call the quantity div(Vw/ |Vm|) the 
nonlocal curvature, and we will denote it A for short. 

Main result. The goal of this paper is to introduce a notion of viscosity solutions of ( |1.1[ ), 



including (13 1, and prove its well-posedness as stated in the following theorem. 



1.1 Theorem (Main theorem). Suppose that a continuous function F : R n xR — » R is degenerate 
elliptic in the sense of (13). Then the initial value problem < \ 1 . 1 | )-( fL2| ) with uq £ C(T") has a 



unique global viscosity solution u £ C(T" x [0, oo)). If additionally u Q £ Lip(T n ), i.e., u is a 
periodic Lipschitz function, then «(•, t) £ Lip(T n )for all t > and 

||V»(-,f)IL<l|VKolL fort>0. 

We also establish a comparison principle for viscosity solutions, and a stability of solutions 
under approximation by regularized problems. Consequently, we observe that in the case of the 
standard total variation flow equation ( |1.4[ ) our viscosity solutions coincide with the semigroup 
(weak) solutions given by the theory of monotone operators. 

Viscosity solutions. The theory of monotone operators does not apply to problems that are 
not of divergence form. Fortunately, the problem ( |1.1[ ) still has a parabolic structure so that 
it is expected to have a comparison principle or an order preserving property. This structure 
puts < | 1 . 1 b in the scope of the theory of viscosity solutions. However, the conventional viscosity 
theory for degenerate parabolic problems does not apply to < | 1 . lb because the quantity 

,V«\ 1 
A = div — — := — trace X\I 



a 



|Vw|/ \q\ \ \q\ 

is unbounded around q = 0, where q = Vu and X = V 2 u, the Hessian of u. In fact, the singularity 
at q = is so strong that the problem has a nonlocal nature in the sense that A is a nonlocal 
quantity on the place where Vu = 0. This phenomenon still occurs even if ( |1.1[ ) is of a divergence 
form p3| so the equation like ( |1.1[ ) is often called a very singular diffusion equation p0| . If the 
unboundedness near Vu = is relatively weak so that the problem is still a local problem like 
jc-Laplace diffusion equation u t - d\v(\Vu\ p ~ 2 Vu) = with 1 < p < 2, the theory of viscosity 
solutions has been well-established [ [3~3~||42l|43||50] |. There has been a considerable effort to 



generalize the viscosity theory to ( |1.1| ); however, so far the results have been restricted to the 
one-dimensional case p4l|2~8}|3~T| or to related level set equations for evolving planar curves 



1 27 28 1; see also review paper [32 1 



We give a new notion of viscosity solutions so that the initial value problem is well-posed in 
higher dimension at least for periodic initial data. Following the previous results, we adjust the 
class of test functions to be compatible with the strong singularity of the operator at Vu = 0. To 
handle this situation, we generalize the notion of a facet, which will include arbitrary closed sets 
with nonempty interior. Then we carefully craft an appropriate class of admissible faceted test 
functions, which can be viewed as a natural extension of the one-dimensional case [|24). This 
class consists of Lipschitz functions whose level set is a facet with a smooth boundary. The 
value of the nonlocal curvature div(V0/ |V0|) for such a (p can be interpreted as the divergence 
of a vector field that is a solution of a certain vector- valued obstacle problem on the facet. In 
contrast to the one-dimensional case, the nonlocal curvature might be non-constant and even 
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discontinuous (see |10l[T3]). Nevertheless, we observe that it is only necessary to evaluate its 



essential infima and suprema on balls of small radius, and a pointwise definition of viscosity 
solutions is possible. We finally finish the discussion of the construction of the class of test 
functions by pointing out that the operator is a standard differential operator when the gradient of 
the solution is nonzero and the conventional test functions of nonzero gradient can be employed. 
The two crucial steps in establishing a reasonable theory of viscosity solutions are proving the 



comparison principle (Theorem 4.1 ) and proving the existence result (Theorem 5.7 ). The proof 
of the comparison principle is quite involved because of the nonlocal nature of the curvature-like 
quantity A = div(V0/ |V0|), which may be discontinuous on facets. When 4>\,(f>2 are C 2 near 
x and V(pi(x) ± 0, the maximum principle asserts that if (pi < (p 2 near x and (f>\{x) = faix), 
then Ai(jc) < A 2 (x) where A,- = div(V0,/ |V0,|). This type of monotonicity is necessary for 
establishing the standard comparison principle for viscosity solutions [33 1. Our key observation 
is that this type of monotonicity of A is still valid on a faceted region (Theorem |2. 15 ). The proof 
of the comparison principle is based on a contradiction argument by assuming the existence 
of sub- and supersolutions u and v, respectively, satisfying u < v at t = but u > v at some 
point for t e (0, T). We double the variables with an extra small parameter £ e W and study the 
maxima of 



u(x, t) - v(y, s) - »P f 
\x-y-£\ 2 /2s + S, 
\t-s\ 2 /2o- + y/(T-t) + y/(T-s), 



for x,yeT\t,s e[0,T], 



where s, cr, y are small positive parameters. The introduction of £ helps flatten the profiles of 
u and v near a maximizer (x, i, x, s) of <Do. The history of this device, which was developed in 
p4|, goes back to JT7| and [|42|. Excluding the case where the conventional theory applies, we 



may conclude (Corollary 4.5 ) that, by fixing small parameters e, cr, y, 



u(x, t) < u(x, t) + S (t, s) - S (i, s) 
v(y, s) > v(x, s) + S (i, s) - S (?, s) 



for x 6 V\ 
forv e U\ 



for t,s e [0, T] for some A > 0, where 

U = {x e T" : u(x, t) > u(x, t)} and V = {y e T" : v(y, s) < v(x, s)}, 

and U A and V A are ^-neighborhoods of U and V, respectively. This flattening procedure gives 
us enough room to construct the desired admissible test functions for both u and v so that their 
related nonlocal curvatures are ordered, which yields a contradiction. 

The existence result is proved via stability under approximation by smooth problems (Theo- 
5.2 ) and by constructing barriers to prevent initial boundary layers (Lemma |5.10 ). The proof 



rem 



of stability is nontrivial because of the discrepancy of test functions of the limit problem and 
of the smooth problems. The key idea is to approximate the spatial profile / of a test function 
using the resolvent operator, that is, formally, 



f a = (I + adET l f fora>0, 



where E is the total variation energy and dE is its subdifferential. In our procedure, we 
approximate E by a smooth energy E m to approximate f a by 



Sam = i 1 + adE m ) f. 
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The merit of the resolvent operator is that the nonlocal curvature is well approximated by 
(fa,m - f)/a, and f a<m is a suitable test function of the smooth problem. The construction of 
barriers is also nonstandard but it is possible to elaborate the idea of |26l. 

We are fully aware of the importance of other boundary conditions like the Dirichlet condition 
or the Neumann condition. However, we do not intend to study them in this paper as these 
problems seem to need several new ideas. 

Literature overview. There is a large amount of literature devoted to the study of the motion 
by crystalline mean curvature. They are classified into at least three approaches-polygonal, 
variational and viscosity. 

The first approach-a polygonal approach-is to understand the crystalline curvature n y for a 
class of crystal bodies like polygonal shapes. For two dimensional crystals, k 7 is a constant on 
facets parallel to the facets of the Wulff shape W T Therefore it is possible to define evolution of 
solutions for a class of polygons having such facets when y is crystalline in strict sense so that its 
Wulff shape is a convex polygon. Such a special family of solutions is called a crystalline flow 
or a crystalline motion introduced in [|2j 52 1. Although there is a higher dimensional approach 
1 35 1, its validity is limited since the anisotropic curvature k 7 might not be constant on flat facets 
1 10 1 and instant facet breaking or bending may occur. For a further development of crystalline 
motion the reader is referred to [(44). 

The second approach-a variational approach-is to understand k 7 as a subdifferential of the 



corresponding singular interfacial energy when the equation has a divergence structure. In [23| 
it is shown that crystalline motion can be interpreted as the evolution given by the abstract theory 
of monotone operators p4}|47j . This approach enables us to approximate crystalline motion by 



evolution by smooth energy and vice versa; see also [39 40] for approximation by crystalline 
motion. However, even for one -dimensional problem it is noted that k 7 + c in ( |1.6[ ) may not 
be constant when c is not spatially constant if the solution is given by the abstract theory [25 1. 
This situation is important since c is often spatially non-constant in the theory of crystal growth. 
In this case if one takes the first approach for planar curve or graph-like function, polygon is 
not enough to describe the phenomena. However, if one allows to include bent polygons with 
free boundaries corresponding to the endpoints of a facet, it is possible to give a rather explicit 
solution (34N37H38). 



The situation becomes much more delicate for higher dimensional crystals even if there is 
no c. As we mentioned before, the anisotropic curvature k 7 might not be constant on facets, 
and instant facet breaking or bending might occur [ [TO") . This suggests that one cannot restrict 
the class of solutions only to bodies with flat facets like polytope with facets parallel to those 
of the Wulff shape in dimension higher than one or in the case c is not spatially constant. A 
more general theory is necessary. A notion of generalized solutions and a comparison principle 
was established through an approximation by reaction-diffusion equations in [[8}|9]| for V v = yK y . 
However, the existence is known only for convex compact initial data [17). 

In general, the curvature k 7 is known to be only of bounded variation and bounded on a 



facet [ 11 , 12 1 . One has to solve a nontrivial obstacle problem to calculate k 7 as noted in [36]. 
The facets where k 7 is constant are called calibrable |T3) . The study of this property goes 
back to a work of E. Giusti [41], where a sufficient condition for calibrability is given for a 
two-dimensional convex domain. Calibrability is related to the study of what is called Cheeger 



sets [45 1 . 



The third approach is an approach based on the theory of viscosity solutions that relies on a 
comparison principle structure of the problem which we already discussed. This is the approach 
taken in this paper. The merit is that one can prove existence and uniqueness in the same class 
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of solutions for ( |1.1[ )-( [L2] ). The viscosity approach and the variational approach for equations 



of divergence form are compared in a review paper [29|. 

There is another approach to definition of solutions of via an original definition of 
composition of multivalued operators which, while restricted to one dimension at the moment, 
allows for the study of facet evolution for quite general data as well as the regularity of solutions, 



including the forcing depending on x and t [46 49 1 



Outline of the paper. We begin the exposition in Section [2] where we give a rigorous inter- 
pretation of the nonlocal curvature operator div(Vw/ |Vw|). This development is then used in 
Section[3] in which we introduce the new notion of viscosity solution. We proceed with the proof 
of the comparison principle for viscosity solution in Section |4} Finally, Section [5] is devoted to 
the proof of existence of a unique viscosity solution of the problem via an approximation by 
regularized parabolic problems. 



2 Nonlocal curvature 

In this section we will discuss the interpretation of the term div (Vw/ |Vw|) as the subdifferential 
of the total variation energy functional. 

2.1 Total variation energy 

The interpretation of the term div (Vw/ |Vw|) is motivated by considering the total variation flow 
in a smooth bounded domain Q. in R" with the Neumann boundary value condition. Its formal 
form is 



w ; = div(^,) in nx(0,oo), 



|=0 on dftx(0,oo), (2.1) 

w| f=0 = uq on Q, 

where v = v^n denotes the unit exterior normal of dQ. so that |^ is the exterior normal derivative 
of u. Since the equation p.l) is singular at Vw = 0, it is nontrivial to formulate a solution 
in a rigorous way. Fortunately, the classical theory of monotone operators applied to the 
subdifferential of a total energy gives a rigorous interpretation [1]. Let us briefly recall its theory. 
We consider the total variation energy E = E n , 



E(ifr) 



f/jV^I if e BV(Q) n L 2 (Q), 
oo otherwise, 



in the Hilbert space H = L 2 (Q). The energy E is convex and lower semi-continuous on H 
and evidently £ J oo. Hence a general theory of monotone operators, initiated by Y. Komura 



[47 1 and developed by H. Brezis [14| and others, guarantees that there is a unique solution 
if/ : [0, oo) — » H which is continuous up to t = and absolutely continuous in [S, T] for any 
T > 5 > such that 

du 

— 6 -dE{u{-, 0) fora.e. t > 0, and w|, =0 = w (2.2) 
at 

for a given arbitrary uq e H. Here dE denotes the subdifferential of E, i.e., 
dE(iff) = {w e H | E(i/f + h) - E(fff) > (w, h) for all h e H}, 
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where < , ) is the inner product of H. In our setting it is the L 2 -inner product (w, h) = f Q wh dx. 
It can be shown that dE(if/) is a convex, closed set and in general it can be empty depending on 
if/ e H. We introduce the domain of dE 

D(dE) = {if/<=H: dE(if/) ± 0}. 



The equation (2.2) is a rigorous form of (2.1 ). Note that the solution u is unique although the 



equation (2.2 ) may seem ambiguous. It turns out that the solution u is right differentiable for all 
t > and its right derivative d + u/dt agrees with the negative of the minimal section {canonical 
restriction) d°E(u(t)) of dE(u(t)). In other words, 



d + u 
~~d7 



-d°E(u(; 0) for all t > 0, 



where 



d Q E(if/) := argmin j|M|# | w e dE(if/)}. 



Since dE(if/) is closed convex, d Q E(if/) is uniquely determined. There is a characterization of 
the subdifferentials which can be found in [[TJ Lemma 2.4] or [|4j. It states that w e dE(if/) for 
if/ e L 2 (Q) if and only if if/ e B V(Q.) and there exists z e L°°(Q.; W) such that 



llzll 
z ■ v m 



w = -divz e L 2 (Q), 
< 1 and (w, if/) = E(if/), 
on dQ (in H~ 1/2 sense), 



(2.3a) 
(2.3b) 
(2.3c) 



where v^n is the unit outer normal vector to dQ. In particular, if if/ e Lip(Q) then Vif/ e L°°(Q) 



and the condition ( |2.3b[ ) is equivalent to z(x) e dW(Vif/(x)) for a.e. x e Q., where dW is the 
subdifferential of W(p) = \p\. It is given as 

AWf , pi(0) P = o, 



ll/'l 



} p^o, 



where 5 r (jc) denotes the open ball of radius r centered at x We observe z(x) = Vif/(x)f \Vif/(x)\ 
whenever Vif/(x) ± 0. 

To motivate the following construction, let us heuristically study the canonical restriction 
d°E(if/) for a simple profile. Assume that if/ takes its maximum in a closure Q of a dom ain 
with the property that Q_cD and if/ < in Q\Q.-. Assume that z satisfying ( 2.3a )-( 2.3c I exists 
and z = Vif/(x)/\Vif/(x)\ in Since if/ attains its maximum on z ■ v^n_ = -1 on at 

least formally, where is the unit exterior normal of dO. , and we observe that 



-d°E(if/) 



div zq 



I div 



IVtAI 



in Q_, 
in Q \ 



such that div zq minimizes the variational problem 



mm 



I div z\ dx 



z ■ v dn _ 



= -1 on dQ, and |z| < 1 a.e. in Q_ 



(2.4) 



The quantity div (Vif/f \Vif/\) depends on the second derivative of if/ and it is nonpositive if if/ is 
concave. This quantity is formally the one-Laplacian of iff. Since it is not a local quantity on 
Q_, we rather would like to call divz a nonlocal one-Laplacian of if/ on Q, . Geometrically 
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speaking, it corresponds to a nonlocal mean-curvature-like quantity. In this paper we shall call 



this quantity a nonlocal curvature for short. The variational problem (2.4) is still convex as it 
admits at most one minimizer div z, although there are several degrees of freedom in the choice 
of zq. Note that divz does not depend on the choice Q. provided that Q c Q. In the literature 
the vector field z satisfying z • v d[i _ = - 1 on dQ.^ and |z| < 1 is often called a Cahn-Hoffman 
vector field. 

For our purpose, we have to consider the total variation energy for a periodic function. 
A periodic function in R" is a function satisfying (p(x + e ; ) = tff(x) for all x e R" and some 
orthogonal basis {ek}'l =i of R" (independent of x). It is identified with a function on the n- 
dimensional torus T" = R'/T where Y is the lattice Zei + • • • + 2e„. The set T" is interpreted as 
the set of equivalence classes {x + Y : x e R"} with the induced metric, topology, and absolute 
value, namely 

dist(;c,y) = dist Rn Qc + T,y + Y), \x\ : = dist(y, 0) = inf \x + £| RB , x,y e T. (2.5) 

In particular, a ball B r (x) is defined as B r (x) := {y e T" : \x - y\ < r). To simplify the exposition, 
we set r = Z" throughout the paper. 

We consider the total variation energy E = E v , for periodic L 2 functions. It is of the form 

F , n ffclViAl if^eBV(T")nL 2 (T n ), 

Ejn(tff) = < J!L 

oo otherwise, 



in the Hilbert space H = L (T). The unique solvability of ( |2.2| ) as well as the characterization of 



the speed is still valid for E = Ej« with obvious modification. For example, ( |2.3c[ ) is unnecessary. 



The characterization ( |2.4[ ) remains unchanged provided that Q. is contained in a fundamental 
domain of T". In this case, -d°Ejn(i//) agrees with -d°E(ifr) in £1 . In what follows, we always 
take E = E v instead of E^. 



2.2 Resolvent equation 

It is possible to approximate the evolution of ( |2.1[ ) via an implicit Euler scheme with a time step 
a > that has the form of the following resolvent problem on T": for given a > and / e L 2 (T") 
find u e L 2 (T") such that 



u + adE(u) 3 f. 



(2.6) 



Due to the standard theory of monotone operators, this problem has a unique solution 
u 6 L 2 (T"), i.e., for every / 6 L 2 (T") there exist unique u e D(dE) and v 6 dE(u) such that 
u + av = f; see for instance [21 §9.6.1, Theorem 1]. 



The following comparison theorem was proved in [ 16 1. 



2.1 Theorem. Let U\,u 2 e L 2 (T") be two solutions of ( |2.6[ ) with right-hand sides fi,fz 6 L°°(T"), 
respectively. If f\< f%, then u\ <u%. 

Proof. We realize that the solution u e L 2 (T") of ( |2.6[ ) provided by the theory of monotone 
operators is a solution in the sense of [[16, Definition 2], and therefore |[T6j Theorem 2] applies. 

□ 



Because dE is a maximal monotone operator, the resolvent problem ( |2.6[ ) offers an alternative 
to finding the minimal section element of the subdifferential dE. 
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2.2 Proposition. Whenever f e D(dE), we have 

■ f " ~ 1 ■ -d°E(f) in L 2 (T) asa^O, 



a 

where f a is the solution u of ( |2.6| ) for given a > 0. 



Proof See [0 Theorem 3.56] (also [(21} §9.6.1, Theorem 2]). □ 



2.3 Nonlocal mean curvature of a facet 

To properly capture the evolution of "non-convex" facets, we need to investigate a more general 



profile than the one considered in (2.4). A pair of open sets (f2_, Q + ) is all that is ne eded to 



describe a general facet D = Q_ \ Q + . Indeed, it will be shown below (Proposition 2.9 ) that the 
quantity div(Vt/r/ |V^|) on the facet D of ft depends only on the sign of ft outside of the facet D. 
This motivates the following two definitions. 

2.3 Definition. Let c T" be a pair of open sets such that Q + c Q_ and Q_ \ Q + ± 
0. We say that (Q_,Q + ) is a smooth pair if Q+ have C 2 (possibly empty) boundaries and 
dist((9Q_, 3Q.+) > 0. This includes the cases when £1+ = or Q- = T", i.e., when one or both of 
dCl+ is empty and the distance from an empty set is +oo by definition. 

2.4 Definition. Let Q+ be a pair of open sets, D.+ c Q- c T", and ft e Lip(T"). We say that ft is 

a support function of(Q-,Q+) if 



ft(x) 



' > in Q + , 

= on D :=Q-\ Q + , (2.7) 
< in Q_ , 



where Q_ := T" \ Q_ denotes the complement of the closure of Q . 

2.5 Remark. It is easy to see that if ft is a support function of a smooth pair (Q.^, Q + ) then -ft is 
a support function of the smooth pair (Q + , Q_). 

For an open set U c T" we define [[3j 

X 2 (U) := [z £ L°°(U;W) : divz 6 L 2 (U)}. 

If dU is Lipschitz continuous, the normal trace z ■ Vdu is well-defined in the sense of H~ l/2 (dQ.) 
for every z 6 X 2 (U), see e.g. [22]. However, here we rather use an equivalent definition based 
on Green's formula. The following version of Green's theorem for vector fields in X 2 (U) was 
proven in 01 for bounded open sets U c R", but the modification for T" is straightforward. 

2.6 Theorem. Let U c T" be an open set with Lipschitz boundary and let z 6 X 2 (U). Then 
Z-Vgu e L°°(dU) such that \\z ■ v du \\ L «, (dU) < \\z\\ L °° m , and 

J wdivz+ J z ■ Va = J u(z ■ v du ) d^H"' 1 for all u e Lip(U). 

JU JU JdU 

2.7 Definition. Let (Q_,D + ) be a smooth pair. We say that a vector field z e X 2 (intD) for 
D = Q_ \ Q + is a Cahn-Hoffman vector field in D if z satisfies 

|z| < 1 a.e. in D (2.8) 
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as well as the boundary condition 

z ■ v aa _ = -1 on dQ.-, z ■ (-v d nj = 1 on dQ. + , 
a.e. with respect to < H n ~ l . 



(2.9) 



The following proposition shows that the boundary condition ( |2.9[ ) is natural and indeed 
holds for any vector field z e X 2 (T n ) that satisfies ( |2.3b| ) for a given ij/ e Lip(T") with a 
sufficiently smooth level-set. 

2.8 Proposition. Suppose that u e Lip(T") and define G t := {x e T" : u(x) > t). Suppose that 
dGo 6 C 1 and that there exists z g X2(T n ) such that z e dW(Vu) a.e. on T". Then z ■ vgG = -1 
a.e. on dGo with respect to < H n ~ x . 



Proof. Since ||z|L < 1, Theorem 2.6 implies that z ■ v dGo > -1 a.e. with respect to IH' 1 and 



f divz= f z-y 5Go d9i n - 1 > -9i"-\dG {) ). 

J Go J dGo 



(2.10) 



Hence it is enough to prove 



f divz<--K" 

-'Go 



l (5G ). 



(2.11) 



We define the cutoffs 
Clearly 



uu = min(max(w, 0), k), k>0. 



( Vu(x) a.e. x, < w(x) < fc, 
1 a.e. otherwise. 



In particular, z & dW(Vu k ) a.e. and therefore 

f z-Vu k = | 



Theorem 2.6 on the left and coarea formula for BV function uu on the right imply 

- f u k divz= f z-Vu k = f \Vu k \= f f |V^ G( | rfxA, (2.12) 
Jt" Jt" Jt" Jo Jt" 

where xg, is the characteristic function of G t . Using Fubini's theorem, we rewrite the first term 

as 



u k (x) div z(x) dx 



i'" 




div z(x) dt dx 

I I div z{x) dx dt. 
Jo Jg, 



IT" JO 



(2.13) 



Since (|2.12|) and (|2.13|) hold for all k > 0, we conclude that 

a.e. t > 0. 



f divz = - f |y^ c 
Jg, Jt" 



(2.14) 
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Clearly G t /* G as t \ and therefore the dominated convergence theorem implies 



J div z = lim I div z, 
Go ^~ * ~*~ *JGt 



and^- G/ — > ;f Go in L^T"). Let f m be a sequence of ? for which ( |2.14| ) holds and t n 
The lower semi-continuity of total variation implies 



as m — > oo. 



f divz=lim f divz = -lim f |V* G J < - f |V^ Go | = - r H n ' l (dGo), 

JGo m ^°° Jg„„ " wo ° Jt» " JGo 



where the last equality holds since dG eC . Thus we have established ( |2.1 1| ) and the lemma 
follows. 

We consider the variational problem 



□ 



mm 



I 



w\ 2 dx 



w 



div z and z is a Cahn-Hoffman vector field in D 



(2.15) 



Of course, (2.4) corresponds to the case when Q + is an empty set. Since (2.8) and (2.9) are 
fulfilled for convex combinations of Cahn-Hoffman vector fields as well, the problem ( |2.15[ ) 
is a convex (obstacle) problem. Moreover, it is lower semi-continuous in L 2 (D) with respect 
to w (the set of Cahn-Hoffman vector fields is closed in L 2 (D)). Thus there exists at least one 
minimizer w provided that there is a Cahn-Hoffman vector field. Fortunately, if the pair (Q_, Q.+) 
representing the facet is smooth, such a vector field exists. Let us recall that w e -dE{ij/) for 
if/ e Lip(T n ) if and only if there exists z e X 2 (T") such that 

divz = w and z £ dW(Viff) a.e. on T". (2.16) 

2.9 Proposition. Let (f2_,f2 + ) be a smooth pair. Then there exists a support function i]/ of 
in the sense of Definition 2.4 such that i]/ e D{dE), and a Cahn-Hoffman vector field 
in D = Q_ \ D. + . Consequently, there exists Wq, the unique solution of the variational problem 
( 2.15 ), andwo(x) = -d°E(i//)(x) a.e. in D. In fact, the value of —d° E(ifj) on D is independent of 
the choice of a support function ip of(£l-, Q + ) such that if/ 6 T)idE). 

Proof. For a given nonempty set E c T", let us define the signed distance function 

d E (x) = dist(jc, E) - dist(x, E c ). 

Since (Q_,Q + ) is a smooth pair, 6 := dist(<9Q_, 3Q.+) > 0. Moreover, the boundaries 
dQ ± are compact C 2 surfaces and therefore they are sets of positive reach: there exist positive 
constants 5-, 6 + , such that da ± e C 2 in {x : dist(x, 8Q.+) < 6+}; see [ 19 48 1. 

We set 5 := | min(5_, 5 + , do) > 0. Let 9 e C°°(R) be a smooth nondecreasing function such 



that 



< 9\a) < 1 for all cr e 

We also define f e Lip(R), 



1 1 

£(cr) := max(min(cr, -8), --6) 



6(o-) 



cr 



for \o-\ < 5/2, 



\8 signer for \cr\ > 8. 



for \cr\ < 8/2, 



^signer otherwise. 
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Finally, we can conclude the construction by defining the support function if/ e Lip(T"), 



-£(d a+ ) in Q+, 

inD :=Q-\a + , 

-£(d a _) in 



and the vector field 



z(x) = -V [9(dnXx)) + 9(dnM)] = -0'(d n _(x))Vd Q _(x) - &(d a+ (x))Vd^(x). 

By construction, z e C l (T n ) and therefore z e X 2 (T"). Indeed, note that 6'(d Q± ^)) * only in 
the (^-neighborhood K+ = [\d a± \ < s\, and d n± e C 2 (K ± ). Moreover, as K + n K. = 0, & < 1 and 

|VJn ± | ^ 1 a.e., we have |z| < 1 in T". This also shows that z satisfies the boundary condition 
< |2.9| ) and thus z is a Cahn-Hoffman field in D. In particular, the variational problem ( |2.15[ ) has a 
unique solution w e L 2 (D). 

Finally, a simple calculation yields that z = Vifj and \z\ = \Vif/\ = 1 whenever Vifr £ 0. Hence 
(TTToT) holds and therefore ^ e D(dE) and divz e -dE(ifr). 



Let wo be the solution of the variational problem ( |2.15[ ). There exists a Cahn-Hoffman vector 
field Zo in D such that w = divz m D. Let ^ be an arbitrary support function of (f2^,Q+) 
such that if/ 6 D(dE). Again, there exists 6 X 2 (T n ) which satisfies ( |2.16[ ) with w = -d°E(if/). 
Clearly z^ e X 2 (intD) and satisfies the boundary condition ( |2.9[ ) due to Proposition [278] and it 
is therefore a Cahn-Hoffman vector field in D. Moreover, the vector field 



z(x) 




on D, 
on D c 



belongs to X 2 (T"), and satisfies the inclusion in ( |2.16[ ). Therefore ||divz^|| L2(£)) < ||div Zol lz. 2 (Z5> 
since divz^. = -d°E(ifr). But the uniqueness of the variational problem yields 

-d°E(ifr) = div z,/, = w a.e. on D. 



This holds for an arbitrary support function if/ e D(dE). 



□ 



2.10 Remark. It is easy to see that if if/ e D(dE) is a support function of a smooth pair (Q.-, Q + ) 
then, for any strictly increasing 6 e C 1 ((-c»,0]) n C l ([0, oo)), Q o if/ e D(dE) with the same 
Cahn-Hoffman vector field zo as if/ on Q._\Q. + , and o if/ is a support function of 12+). This 
fact will be used multiple times with the particular choice 9(s) = as + - fis- for some a,/3 > 0. 

We are now ready to define the nonlocal curvature for a smooth pair (Q_, Q + ). 



2.11 Definition. Let be a smooth pair. Let wq be the unique minimizer of ( |2.15[ ). Then 

we call wo the nonlocal curvature of f2 + ), and we denote it by A(£X_, 

Although there is a huge literature on obstacle problems (see e.g. pTJ), for our obstacle prob- 



lem ( |2.15[ ) there seems to be no literature when D + is non-empty except in the one-dimensional 
case in which wo is always constant. (If one end of f2_ is not covered by f2 + , w must be zero 
near that point.) In general, it is hard to calculate w . If there is no obstacle condition \z\ < 1 
in (|2.8|), the minimizer w is always constant since its Euler-Lagrange equation is Vdivz = 0. 
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If the minimizer w is a constant A, the value of A is easily computed by the Green's theorem, 
Theorem |2 .6[ Indeed, 

<H n (D)A = f div z dx= f z -v dn _d'H n - l + f zo ■ (~v d£1+ ) d'H"- 1 

JD JdSl- JdSl+ 

= - < H n -\dQ-) + 'H n -\dQ. + ) 

and therefore 

A= ^)(' H "" 1 ^- <Hn " 1 ^ Q -))' 

where < H n denotes the rc-dimensional Hausdorff measure so that < H n is the Lebesgue measure. 

The question is whether or not w is constant on D. When Q + is empty, this problem has 
been well-studied by Bellettini et al. [ fT3| . A facet where wo is constant is called calibrable. 
There is a necessary and sufficient condition for calibrability. A sufficient condition for a convex 
set to be calibrable is that the principal curvature of <9Q- is bounded by r H n ~ l (dQ.^)/'H"(D). 

In the case when w is not constant, not much is known except regularity [ |TTj[T2j . Among 
other results, they proved that the minimizer is in 5V(fi_) but may be discontinuous. See also 
© Theorem 5.16]. 

Let us point out here that our (Q_, Q.+) pair formulation does not require that D is connected 
so it is more general. 

We conclude this subsection by illustrating the definitions on a few examples of simple 
smooth pairs. 

2.12 Example. Let Q = T" and Q.+ = 0, i.e., the facet D = Q_ \ Q + is the whole torus. In this 
case A(Q_, Q.+) = on D = V since z = is a Cahn- Hoffman vector field in D. 

2.13 Example. Let Q_ = B R (0), R e (0, 1/2), and let Q.+ = 0. Then 



is a Cahn-Hoffman vector field on D. Note that the minimizer is wq 



—n 



2.14 Example. Let Q+ = B R± (0) with < R + < i?- < 1 /2. The minimizer wq of ( |2.15| ) for this 
pair Q + ) is the constant 

A(H_, Q + ) = (-K"- 1 ^,) - ^"- 1 (5Q_)) . 

In particular, the facet D is calibrable even though f2 + ^ 0. 

Proof. We shall construct a Cahn-Hoffman vector field z so that div z is a constant. In fact, we 
take z of the form z = Vip, where tp is a radially symmetric solution of the Poisson equation with 
Neumann boundary data. Let us denote 

?fl— 1 D'l-l 



is a solution of 



:= A(fl_, Q + ) = n— 



\A<p = b in D, 

lg = ±l onW=/? ± . 
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In fact, (p might be chosen in the form 



a i b n 

<p(x)=~ \x\ 2 ~ n + -\x\ 2 . 

2 - n 2 



To satisfy the boundary condition, we take 

R'!_ - R'l 



A simple computation yields that 

\x\ 

where 6(\x\) = a \xt n + b\x\> 0. Since 6{p) is convex for p > and 6(R±) = 1, we conclude 
that \z\ e (0, 1] in D. Therefore z is a Cahn-Hoffman vector field. □ 



2.4 Monotonicity of nonlocal curvatures 



In this subsection we shall prove the monotonicity of nonlocal curvature A with respect to a 
smooth pair Q+). Its precise form is given by the following theorem. 

2.15 Theorem (Monotonicity). Let Oi + ) be a smooth pair for i = 1,2. Let A' = A(QL, ^ ! +) 
be the nonlocal curvature on D l = Cl' + for i = 1,2. If£l l _ c Q 2 and Q. l + c Q. 2 , then A 1 < A 2 
a.e. in D l n D 2 . 



Proof. We apply the comparison principle for the resolvent problem ( |2.6| ); a similar idea was 
pursued in [35 1, where the evolution equation <\2.2} is used instead of the resolvent problem. 
We have 

A'' = -d°E(iff') a.e. in D\ 

where if/' e D(dE) are the support functions for (Q'_, £l' + ) constructed in Proposition 
that it is possible to take the same 6 > in the proof of Proposition 2.9 for both iff 
then by construction and due to the ordering Q,]. c Q 2 we have 



2.9 Observe 
and if/ 2 , and 



•A 1 <<A 2 



in 



iff 1 = if/ 2 = 



inZ) 1 C\D 2 . 



(2.17) 



For each a > 0, we consider the solution f l u of the resolvent problem ( |2.6[ ) with right-hand 
side iff 1 . Due to Proposition 2.2, (f l a - iff 1 ) / 'a — » -d°E(iff l ) in L 2 as a — > 0. Therefore there is a 
subsequence a* — > as — > oo such that (f l clk - if/ 1 ) I at — > -d°E(if/') a.e. as k — > oo for z = 1,2. 
The comparison principle, Theorem 



D 1 n D 2 for all it. Therefore 



2.1 



and (|2T7J) imply that £ < / 2 and £ - ^ = £ in 



a* 



>a k 



A 1 = -d°E(iff Y ) = lim ^ < lim = -d°E(if/ 2 ) = A 2 a.e. in D 1 n £> 2 

and the monotonicity of A is established. 

To illustrate the previous theorem, we show some basic bounds on A(Q_, Q + ). 
2.16 Proposition. Let (Q_, Q + ) be a smooth pair. Then 



□ 



|A(Q_,Q + )(*)| < 



p(x) 



for a.e. x e D := Q- \ Q.+, 



where p(x) := sup |p < 1/4 : x e 5 p (j) c Dfor some y\. 
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Proof. Let iff be the support function of (Q_, Q + ) from Proposition [Z9] and set M = max iff 



We fix £ e D and choose p e (0, 1/4) and y for which £ 6 B p (y) c B p (y) c int D. There also 
exists (5 6 (0, 1 1 A) such that 5 p+(5 (y) c D. 
We introduce the barrier tp e Lip(T"), 



|x-£l<P, 
M ttP p<|*-£|<p + £, 

M \x - £| > p + 6. 



This function is a support function of the smooth pair (' 
Example 



2.13 



A), where A = T" \ B p {£). As in 
it can be shown that <p e D(dE); that example also shows that A(T",A) = - 



on B p (y). Clearly, Q + c A and therefore the monotonicity result in Theorem 2.15 allows us to 



compare A(Q_, f2 + ) < A(T, A) = - a.e. in BJg). The lower bound is analogous. 



□ 



3 Viscosity solutions 

In this section we finally introduce a new notion of viscosity of ( |1.1| ). As in the previous works, 
it is necessary to carefully chose an appropriate class of test functions. To this end we introduce 
admissible faceted test functions. 

3.1 Definition. Function ip(x, t) = fix) + g(t), where / 6 Lip(T") and g e C l ([t -6,t + 6]) for 
some 8 > and i e (0, T), is called admissible faceted test function at (x,i) e Q := T" x (0, T) 
with pair Q+) if / is a support function of a smooth pair (Q_, Q + ) and x 6 Q. \ f2 + . The 
pair is called the pair associated with <p at (x, t). 

3.2 Definition. We say that an admissible faceted test function (p at (x, t) e Q with pair (Q.,Q. + ) 
is in general position of radius n > wzY/? respect to u : Q —> R at (x,t) if, for all /* 6 5^(0), we 
have x + h e Q \ Q + and 

w(x -h,t)- <p(x, t) < u(x, t) - (p(x, t) for all x e T, t e [t - rj,t + rj]. 

If such 7] > exists, we simply say that (p is in general position with respect to u at (x, i). 

The notion of general position formalizes for admissible test functions the idea that two 



facets stay ordered in the sense of Theorem 2. 15 even when one of them is shifted by a small 



distance in an arbitrary direction. This makes the family of test functions smaller and turns out 
to be the right ingredient to obtain the existence of solutions through stability, without affecting 
the comparison principle. 

3.3 Definition (Viscosity solutions). An upper semi-continuous function u : Q — > R is a 
viscosity subsolution of ( |1.1[ ) if the following holds: 

(i) (faceted test) If (p is an admissible faceted test function such that (p is in general position of 
radius rj with respect to u at a point (x, t) e Q then there exists 6 e (0, 77) such that 



(p t (x, t) + F\0, ess inf A(£}_,£} + ) < 0, 

\ B 6 (x) J 

where (Q._, D+) is the smooth pair associated with tp at (x, t). 
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(ii) (conventional test) If (f e C^)(U) in a neighborhood U c Qof a point (£, f), such that 
u — tp has a local maximum at (x, t) and |V<p| (x, f) ^ 0, then 

f) + F (V<p(x, t), k(V<p(x, t), V 2 cp(x, t))) < 0, 

where V 2 is the Hessian and 

Itracex(/-^ 
\P\ \ \P\ 



k(p,X) = — traceXj/- for pef"\ {0}, X e 5", (3.1) 



so that k(V(p(x, i), V 2 </?(£, t)) = [div(V^/ |V^|)] (x, t). Here S n is the set of n x n-symmetric 
matrices. 

Viscosity supersolution can be defined similarly as a lower semi-continuous function, replac- 
ing maximum by minimum, < by >, and ess inf by ess sup. Furthermore, in (i) ip must be such 



that -<p is in general position of radius 77 with respect to -u (see also Remark 2.5 ). 



Function u is a viscosity solution if it is both a subsolution and supersolution. 



3.4 Remark. Our definition is essentially similar to that of local version in p4| where one- 
dimensional problem is studied for more general equation u t + F{u x , ( W'{u x )) x ) = with convex 
W whose derivative has a discrete set of jumps, not necessarily a singleton. 



4 Comparison principle 

Our next task is to prove the comparison principle on the space-time cylinder Q := T" x (0, T) 
for any fixed T > 0. 

4.1 Theorem (Comparison). Let u and v be respectively a bounded viscosity subsolution and a 
bounded viscosity supersolution of ( |1.1[ ) onQ.Ifu<vatt = then u <v on Q. 



We will start with a variation of the standard doubling-of-variables argument. As in |24J , we 
define 

w(x, t, y, s) = u(x, t) - v(y, s) 
and, given positive constants s, o~, y and a vector (eT", the functions 

^V ( (x,t,y, s;s,cr,y) := + S(t, s; cr, y), 

2s 

S (t, s; cr, y) := — + - — - + 



2o- T-t T-s' 



where | jc - y - £\ is interpreted as in (2.5 ). 



We analyze the maxima of functions 

® s (x,t,y, s;s,cr,y) := w(x,t,y, s) - Y ? (x,?,y, s;s,cr,y). 

The purpose of introducing an extra parameter £ is to restrict the class of possible contact 
points that require a construction of an admissible faceted test function in general position 



(Definition 3.31)) down to the profiles that allow for small shifts. Indeed, if at least one 
maximum of for small £ occurs at a point (jc, t, y, s) such that x - y £ £ then V^O^ ^ at 
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the point of maximum, and we can apply the classical viscosity argument and proceed with 
the well-known construction of a test function for Definition |33fii). On the other hand, if all 
maxima of (D^ for small £ occur at points with x - y = g, it is necessary to proceed with a more 
involved construction of an admissible faceted test function in general position. Fortunately, 
the structure of maxima with respect to £ provides enough information to infer "flatness" of 
u and v at the contact point, which in turn guarantees that the construction is possible. We 
carefully analyze the profiles of u and v at the level of the contact, construct two smooth pairs 
that represent facets ordered in the sense of Theorem 2.15 and construct admissible faceted test 
function with the help of the support function introduced in Proposition 2.9 The definition of 
viscosity solution then yields a contradiction. 

Therefore, suppose that the comparison principle does not hold, i.e., 

m := sup [u(x, t) - v(x, t)] > 0. 

We have the following proposition: 

4.2 Proposition. Set M := maxg x g w > m , m' Q := |m , k(s) := \ (s(m - m'^j 1 . 

There exist positive constants £o,0"o and y sucn that for any s e (0, so), cr e (0, cr ), 
y e (0, 70) an d £ e T" such that \£\ < k(s) the following is true: 

If (X{, t;,y s , s ? ) e arg maxg x g O f (• ; s, cr, y) then 

(0 (x e ,t e ,y ( ,se) eQxQ, 
(U) \t ( - s c \ < (Mo-) 112 , \x ( -y ( - (\ < (Ms) 1 ' 2 , 
(Hi) sup > m' . 



Proof. See [24 Proposition 7.1, Remark 7.2]. In our case, we take y := min(yo,7^) for 
simplicity. □ 

The possible scenarios can be divided into two cases that will be treated independently: 

• Case I: all the points of maximum lie on the space diagonals, i.e., 

argmax O f c {(x, t,y, s) eQxQ : x - y = ^} (4.1) 

2x2 

for all e, cr, y and |£| < k(s). 

• Case H: there exists e, cr, y, \g\ < k(e) and (x, t,y, s) e arg maxg x g O^- such that x - y ± g. 



4.1 Case I 



Let us fix s e (0, s ), cr e (0, cr ), y e (0, y ) such that ( gTT) holds for all ^, \£\ < k(s). To 
simplify the notation, we set A := k(s)/2. We also do not explicitly state the dependence of the 
following formulas on the fixed s, cr and y. Let us define 



sup <E>^. 

2x2 



The first step is the application of the constancy lemma |24[ Lemma 7.5], as in the proof of 
Proposition 7.6], to show that €(£) is constant for |£| < A. 
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4.3 Lemma (Constancy lemma). Let K be a compact set in R N for some N > 1 and let hbe a 
real-valued upper semi-continuous function on K. Let cpbe a C 2 function on R d with 1 < d < N. 
Let G be a bounded domain in W 1 . For each £ e G assume that there is a maximizer (r^p^) € K 
of 

H ( (r,p) = h(r,p)-<f>(r-0 
over K such that V0(r f - = 0. Then, 

h^) = sup {H ( (r,p) : (r,p) e k) 

is constant on G. 



We apply Lemma 4.3 with the following parameters: 

N = 2n + 2, d = n, p = (y, t, s) e T" X R x R, 
K = {(x- y,y, t, s) : (x,y) e T" x T", (t, s) e [0, T] x [0, T]}, 

G = B 2A (0), 

\r\ 2 

h(r,p) = w(r + y,t,y,s)- S(t,s), (f>{f) = —. 

2s 

K can be treated as a compact subset of W in a straightforward way. We infer that = h^g) 
is constant for |£| < A. 

Therefore we have also an ordering analogous to J24| Corollary 7.9]: 



4.4 Lemma. Let (x, i, x, s) e arg max O . Then 

u(x, t) - v(y, s) - S (t, s) < u(x, t) - v(x, s) - S it, s) 

for all s,t 6 (0, T) and x,y e T" such that \x - y\ < A := k(s)/2. 

Proof For \x - y\ < A, a straightforward computation using the constancy of £ yields 

u(x, t) - v(y, s) -S(t, s) = O v _ v (x, t,y, s) 

< £( X -y) = £(0) 

= u(x,t) -v(x,s) - S(t,s). □ 



We will use Lemma 4.4 to create room for ordered smooth facets, whose support functions 
will play the role of admissible test functions for u and v. The test functions themselves do not 
have to be ordered due to the positive 0-homogeneity of the nonlocal curvature in our setting 
(recall Remark [2 .10| ). 

At this moment, we will fix (x, i, x, s) 6 arg maxg x g O and we set 

a := u(x, t), J3 := v(x, s). 
Let us define the closed sets 

U :={u(;i)>a}, v.= {v(;s)<J3}, 

see Figure [T] 
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Figure 1: Facet construction for u and v in the proof of the comparison theorem. 



It is convenient to introduce the notion of dilation and erosion of a set in morphology. For 
any set A c T" and p 6 R, A p denotes the generalized neighborhood 



A + B p (0) p > 0, 

A _ p = 0, (4.2) 

{jc e A : B|p|(x) c a} p < 0. 



Here A + 5 = {jc + y : x e A, y e B] is the vector sum called the Minkowski sum. In image 
analysis it is often written as A p = A © B p (0) for p > and A p = AG B^O) for p < 0, where © 
denotes the Minkowski addition and © denotes the Minkowski decomposition. In morphology, 
ffi is called dilation and © is called erosion. If A is open and p > 0, we have 

A p = {xeT : dist(x,A) <5}, 
A~ p = {x e A : dist<XdA) > 6}. 

Clearly, for any pel, the generalized neighborhood A p is open (resp. closed) provided that A 
itself is open (resp. closed). Moreover 

A' p = ((A c rr, 

where A c denotes the complement of A, i.e., A c := T" \ A. If A has the interior ball property with 
radius r > then (A' p f = A for any < p < r. 

4.5 Corollary. We have 

u(x, i)<a + S(t,s)-S (t, s) forxeV A , te (0, T), 

and 

v(y, s)>/3 + S (?, s) - S (t, s) foryeU A , se (0, T). 

Proof. If x e V A then there exists y e V with \x - y\ < A and the conclusion follows from 
Lemma K4 and the definition of V. The proof of the second inequality is analogous. □ 

In what follows, we will need to construct smooth sets that approximate a given general set. 
The next lemma provides such approximation. 
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4.6 Lemma. For any set £cT" and constants r > 0, 6 > 0, there exists an open set G c T" 
with a smooth boundary such that 

E r cGcGc int(£ r+5 ). 

Proof. This is a simple consequence of the Sard's theorem. First, if E r+S = T" (resp. E = 0) then 
we choose G = T" (resp. G = 0) and we are done. Thus we may assume that ^ £ c £ , '" +<5 c T 1 . 
Let := dist(x, £) and let 5 /4 be the standard smooth mollifier with support in B S / 4 . Let 
g = d * 4>s/a- We observe that g e C 00 , and, by triangle inequality, also g < r + |# on E and 
g > r + IS on (mt(E r+s )) c . 

Sard's theorem yields that \g({x : Vg(x) = 0})| = 0, where |-| stands for the one-dimensional 
Lebesgue measure. Therefore, there exists p 6 (r + ^6, r + |#) such that Vg ^ on {g = p}. We 
set G = {x : g(x) < p} and observe that dG = {g = p] e C°° by the implicit function theorem. □ 

We continue with the construction of smooth pairs for u and v, both containing x. We set 
r := and consider the closed sets 



Z:=U\V 15r , W:=V\U l5r . 



By Lemma 4.6 there are open sets G„ / = w, v, with smooth boundary such that 

Z 3r c G u cc int(Z 4 '') c Z 4r c U 5r c //„ cc int([/ 6r ), (4.3a) 

W 3r c G v cc int(W 4r ) c W 4r c V 5r c H v cc int(V 6r ). (4.3b) 

If Z (resp. W) is empty then we take G„ = (resp. //„ = 0). 

Because, by construction, Z 4r n y 6r = as well as W 4r n t/ 6r = 0, we also have 

dist(G„, H v ) > 0, dist(// H , G v ) > 0. (4.4) 

Our next task is to construct admissible faceted test functions around (x, t) and (x, s) for u 
and v, respectively, with facets given by the smooth pairs 

(Q_ „, Q + „) := (H u , G u ) and (Q_, v , Q + , v ) := (int G c v , int H c v ). (4.5) 



Note that these smooth pairs satisfy the assumptions of Theorem 2.15 in particular, the facets 
are ordered due to ( |4.4| ). 
Consider the closed sets 



X = {x : distQt, £/) > r}, F = [x : dist(x, V) > r). 

Since w(-, f) < a in X and v(-, s) > in J 7 , and w and v are semicontinuous, there exists r > 
such that 

«(-, < or + ^) - S (t, s) in X for ? £ [f — t, f + t] , (4.6a) 
v(-,s) >/3 + S(t,s)-S(t,s) in Y for se[s-Tj + r]. (4.6b) 

We note that if both X and Y are empty we choose some suitable small r > 0. Therefore, together 
with Corollary |4 .5 [ 

M(-,0<a + S(U)-S(U) in X U V A for t e [t-r,t + T], (4.7a) 
v(-,j)>£ + S(?,j)-S(f,j) inYU f/^forse [j-t,^ + t]. (4.7b) 
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4.7 Lemma. The choice of the sets G[ and H t guarantees that 



(G u y 2r uxuv A = (H u y r ux = r, 

(G v )" 2r U Y U U A = (H v )~ r U Y = T. 



(4.8a) 
(4.8b) 



Proof. We will only show (j478aj), the other equality is analogous. Suppose that there exists 
x £ {G u )~ 2r U X. We will show x e V A . From the definition of X and from ( |4.3a[ ), it follows that 
Z r c (G u )~ 2r and thus 



dist(jc, U) < r, 



dist(jc,Z) = dist(jc, U \ V l5r ) > r. 



(4.9) 



Therefore there exists y e U with \x — y\ < r, and the second inequality in <\4.9) implies that 
y 6 V 15r . The triangle inequality yields 

dist(jc, V) < \x - y\ + dist(y, V) < 16r < 20r = ^. 



The other equality in (4.8a ) follows from (4.3a ), which yields (H u y r UXd U U X = T". □ 



Our goal is to construct admissible faceted test functions that correspond to facets (Q- fi+,0> 
/ = u, v, which are in general position with respect to u and v at (x, t). Therefore, let iff u , iff v be the 



support functions constructed in Proposition 2.9 for smooth pairs (Q_ H , Q.+, u ) and (£l~ iV , n+, v ), 
respectively. We define 

<p u (x, t) = A u [if/ u (x)] + - n u [if/ u (x)]_ + a + S(t, s) - S(t, s), 
<p v (x, s) = A v [i/f v (x)] + - jJL v [if/ v (x)]_ +0 + S(t, s) -S(t, s), 

where [q]+ := max(+<5r, 0). 

4.8 Lemma. There are positive constants A u , /x v , chosen large enough, and positive constants 
H u , Ay, chosen small enough, so that 



u(x, t) < (f u (x - h, t) 



for xeT, te [t - r, t + t], h e B r (0), 



and 



<p v (x -h,t)< v(x, t) 



for xeT, te[s-T,s + r], he fl,.(0). 



Proof. We shall show the statement in the case of <p u . 

Let us define ip r (x, t) = min^,. tp u (x-h, t). It is easy to see that (p r e C(Q) and <//(•, t)-g{t) is 
a support function of the (not necessarily smooth) pair ((H v )~ r , (G u )~ r ). Moreover, the statement 
of the lemma is clearly equivalent to u < (p r for x 6 T, t e [t - r,t + r]. 



The space T" can be divided into three parts with the help of observation ( |4.8a[ ): 



(a) Mi = (X U V A ) C . Since Mi c (G u y 2r c (G u y by ( |48a| ), \j/ u > in G u and (G u y r is open, 
continuity yields min\ h \< r min^- +ft iff u > and we can choose A u large enough to satisfy 

u < ip r in Mi x [i - r, i + r] due to boundedness of u. 

(b) M 2 = (H u y r \ Mi = (H u y n(XU V A ). Observe that ifj u > on H u by construction. 



Therefore u < tp r on M 2 x [t - r,t + t] due to (4.7a). 
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(c) M 3 = ((H u ) r ) c . We can choose ji u small enough so that 



-/•*« ll<AJL ^ max u(x,t)-a-S(t,s) + S(i,s) 

XX[t-T,t+T] 



because the right-hand side is negative thanks to ( |4.6a[ ). Since M 3 c X by ( |4.8a[ ), we see 
that u < (p r with this choice of jx u on M 3 x [f - t, ? + r]. 



The result for <p v can be obtained similarly. 



□ 



Let us recall the definition of smooth pairs (fl_ in ( |4.5[ ). Now, it is clearly true, by 
definition, that ip u (x, t) = u(x, t) and (f v (x, s) = v(x, s). Next, we give a simple argument to show 
that B r (x) c H u \ G u and B r (x) c H v \ G v . This is immediate, by recalling the obvious fact 



x e U n V and its consequence B r {%) c U r P\ V c H u n H v , which, together with ( |4.4| ), yield 
the claim. Consequently, by construction, for all \z\ < r, tp u is an admissible faceted test function 
at (x + z, t) with smooth pair (Q_ M , Q. +:U ) and (p v is an admissible faceted test function at (x + z, t) 
with smooth pair (Q_ v , ^ +jV )- 

Therefore Lemma [4~8] yields that <p u is in general position of radius r with respect to u at 
{x, t), and -(p v is in general position of radius r with respect to -v at (x, 5). 

By definition of viscosity solutions, we must have for some 6 6 (0, r) 

S t (t, s) + F(0, ess inf A(Q_ u , Q + „)) < 

B 5 (x) 

and 

-5 4 (f, s) + F(0,ess sup A(Q_ jV , Q +)V )) > 0. 

B 6 {x) 

Subtracting these two inequalities, we get 

< T^Ti + ,t 7 ^2 ^ i7 (°' ess supA(0_ >v ,a +;V )) - F(0,ess inf A(Q_,„, Q +M )). 



Since Q ± u c Cl ±tV , the facets are ordered and the comparison principle Theorem 2.15 implies 



ess inf A(Q_ „, Q+^ u ) < ess sup A(£l_ v , Q +iV ), 
and we get a contradiction due to the degenerate ellipticity of F. 



4.2 Case II 

Since this case is quite classical, similar to J24| 7.§G], we just give a short outline of the proof. 

Let us first recall that there exist e, cr, y, |^| < and (z, z') e arg maxg x g O f such that 
x-y±£. 



Using the special structure of O^, we can proceed as in the proof of [50 Theorem 3.1], 



using the elliptic result [18, Theorem 3.2] that is then extended to the parabolic result in 
J50| Corollary 3.6], and we obtain symmetric matrices X, Y 6 S n , X < Y, such that 



(^ t ,q>^ x ,X)eP^(u(x,t)), 



(-^, s ,-^, y ,Y)eP^-(v(y,s)). 



Here !P 2 ± are the closures of the second order parabolic semijets, as defined in [18 §8]. 
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Let us observe that 



x — y — £ 
P ■= <&(,x = -®(, y = * 0, 



where p eR" has to be interpreted as the element of the equivalency class x - y - £ + TP with 
the smallest norm. 

Therefore, the definition of the viscosity solution implies 

&!, t + F(p,k(p,X))<0, 
-&^ + F(p,k(p, Y))>0. 

Summing the inequalities, we arrive at 

®(,t + ®C,s ^ F(P, Kp, Y)) - F(p, k(p, X)). 



The right-hand side is nonpositive due to the ellipticity of F, X < Y and the formula (J3TTJ). We 
have 

a contradiction. 



The proof of the comparison principle, Theorem 4. 1 , is now complete. □ 



5 Existence of solutions by stability 

We shall prove the existence of solutions of ( |1.1[ ) via an approximation by parabolic problems 

(u t + F(Vu,-d°E m (u(;t))) = 0, 

["l/=0 = u o, 

where E m are some smooth energies approximating E. For given m > 0, we choose the energy 
E m as a smooth approximation of the total variation energy, given as 



f/ W m (Vifr) dx for if/ e H l (T n ), 
|+oo for e L 2 (T") \ H l (T n ), 

where W m (p) := (\p\ 2 + m" 2 ) 1/2 + m~ l \p\ 2 . 



The main features of the approximation are summarized in the following proposition. Let us 
point out that other approximations are possible, as long as W m form a decreasing sequence of 
smooth, uniformly convex, radially symmetric functions converging to W(p) = \p\ pointwise as 
m — » oo. 

5.1 Proposition, (a) E m form a decreasing sequence of proper convex lower semi-continuous 
junctionals and E = (inf m E m ) t , the lower semi-continuous envelope o/inf m E m in L 2 (T"). 

(b) The subdijferential dE m is a singleton for all i// 6 D(dE m ) = H 2 (T n ) and its canonical 
restriction can be expressed as 

V(/r 



-d°E m (\fr) = div [(V p W m )(V^r)] = div 



+ m l Aif/ a.e. (5.2) 



(c) Due to the ellipticity of F, if —dE m (\ji) is interpreted as \5.2) , the problem ( |5.1[ ) is a 
degenerate parabolic problem that has a unique global viscosity solution for given continuous 
initial data u 6 C(T"). 
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Proof. Clearly W m is a sequence of decreasing functions, and hence E m is a decreasing sequence 
of functionals. The proof of lower semi-continuity and convexity of E m , together with the 
characterization of the subdifferential in (b), can be found in |2TJ. 

To prove (inf m E m ) t = E, let us denote F = inf m E m . We need to show that E{\)j) = F^fa) = 
hm'mf^ k ^ F(fa) for all if/ e L 2 . First, let us assume that \p e H 1 ^"). Then E m (fa) < oo; 
moreover, since W m (D\p) J, W(Dif/) a.e., the monotone convergence theorem implies E(tf/) = 
inf E m (if/) = F{ij/). On the other hand, if iff t H l then E m (fa) = oo for all m and hence F(i//) = oo. 
We therefore conclude that E < F, and, recalling that E is l.s.c. in L 2 , we also see that E < F*. 
Let now iff e BV(T"). There exists a sequence fa e C°°(T n ) c H [ (T 1 ) such that fa -> tfr in L 2 
and F(^fc) = F(^) — > F(i/r); see |[TJ Theorem B.3]. Therefore 

E(ifr) < F*(ifr) < liminf F(fa) = lim inf F(^ fe ) = E(iff), 

k k 

which concludes the proof of (a). 

For (c) see (18). □ 



In the following two subsections we prove that the limit of solutions of (5.1 ) is the unique 



viscosity solution of with the same initial data. First, we establish a key stability property 



of subsolutions and supersolutions of (5.1 ), which we then follow by a construction of barriers 



near the initial data that prevent the formation of an initial boundary layer in the limit. 

5.1 Stability 

Let us recall the definition of the half -relaxed limits 

★-limsup i/ m (jc, t) := lim sup sup u m {y, s), 

m-^oo m >k | y _ v |< l 

\s-t\4 

★-liminf u m (x, t) := lim inf inf u m (y, s). 

m— »oo k— >oo m>k \y-x\<j 



\s-t\<i 



We have the following stability property. 



5.2 Theorem (Stability). Let u m be a sequence of subsolutions of ( |5.1[ ) on T" x [0, oo), and let 
u = ★-limsup m ^ co u m . Assume that u < +oo in T' 1 x [0, oo). Then u is a subsolution of ( |1.1| ). 

Similarly, u = ★-liminf m ^ co u m is a supersolution of ( |1.1[ ) provided that u m is a sequence of 
supersolutions of \5.\\ and u > -oo. 



Before proceeding with the proof of the stability theorem, let us first give a motivation and a 
brief outline. Intuitively, the main challenge in the proof of stability lies in the different natures 



of the nonlocal limit problem ( |1.1[ ) and the local approximate problems ( |5.1[ ). 

The basic idea of the standard stability proof in the theory of viscosity solutions is the 
observation that if u - (p has a strict local maximum at a point (x, t), then u m - <p will have a 
local maximum at a nearby point (x m , t m ), for a subsequence of m. Using the knowledge that u m 
is a subsolution of the approximate problem, classical maximum principle arguments provide 
information on (p at (x m , t m ), and the proof is concluded by taking the limit m — > oo. 

However, the local nature of the approximate problem allows us to recover only local 
information about the test function (p at (x, t). Therefore, in the current setting, this direct 
argument works only as long as V(p(x, t) ± 0. 
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When <p is an admissible faceted test function at (x, t), in particular V(p = in the neighbor- 
hood of (x, t), the viscosity condition on u is formulated via the inherently nonlocal quantity 
A(Q_, Q.+)(x). Fortunately, it is possible to encode the missing nonlocal information into the 
test function itself by means of a small perturbation of tp depending on m (and other extra 
parameters), as long as the perturbation can be shown to vanish uniformly when taking the 
various limits. The perturbed test function method was pioneered in the theory of viscosity 
solutions by Evans p0| , where it was applied to the homogenization of an elliptic problem and 
thus the choice of perturbation was completely different. 

The perturbation (p a , m (x, t) = f a , m (x) + g(0 of ip(x, t) = f(x) + g(t) is motivated by the 
interpretation of div(V// |V/|) as the L 2 -gradient of the total variation energy in the gradient 



flow (2.2) with initial data /. The evolution of the gradient flow can be approximated by the 



implicit Euler scheme with finite time- step a > 0, whose first time-step is in fact the solution f a 



of the resolvent problem ( |2.6[ ). The resolvent problem is, however, conveniently approximated 
by the solution / a m of the resolvent problem for energy E m . Moreover, the nonlocal nature of 
resolvent problems provides the missing information about the shape of the facet of tp. 

We collect the relevant features of the approximation by resolvent problems in Proposition |5.3| 
below. The statement contains an extra parameter e which will be used throughout the proof of 
stability to rigorously justify that the limit problem does not depend on V(p outside of the facet. 

The proof of stability for subsolutions proceeds roughly in the following manner: 1) as test 
functions with Vtp ± can be handled by the classical stability proof, we restrict our attention to 
admissible faceted test functions tp(x, t) = f(x) + g(t) such that u - <p has a maximum at (x, t) on 
a facet; 2) we show that it is possible to assume that the set of points of maxima can be restricted 
to a small neighborhood of the facet; this is a consequence of <p being in general position with 
respect to u; 3) we perturb the test function by solving the resolvent problems for energy E m with 
"time-step" a and right-hand side ef; 4) after using the classical maximal principle arguments to 
infer information about the action of operator ( |5.2| ) on (p am , we sent m — > oo, which provides 
information about the action of the operator ( |1.1[ ) on (p a ; finally, 5) sending a — » and e — > 
together with a geometric argument concludes the proof. 

5.3 Proposition. For f e Lip(T n ) and positive constants m, a, s, the resolvent problems 

f a + adE(fa) 3 ef, (5.3a) 
fas* + adE m {f ajn ) 3 ef, (5.3b) 

admit unique Lipschitz continuous solutions f a and f a<m , respectively, and 

IIV/JL.IIV/^IL^ellV/L. (5.4) 

Moreover, f a<m e C 2 ' a (T n ) for some a > 0. 

In what follows, let e > be fixed and we define 

h ■— £ f h •— f a,m ~ S f — ffij? ( f \ (K Z.\ 

"a • — j n am . — — O tLm\Ja,m)- V -'■-'/ 

a a 

Then, for fixed a > 0, 

f a m =| f a uniformly as m — » oo, and 

h a m =| h a uniformly as m — » oo. 



Additionally, 



f a =^ ef uniformly as a — > 0. 
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If furthermore f 6 !D(dE), then also 



h a -> -d°E(sf) in L 2 (T") as a -» 0. 

Proof. The existence and uniqueness of / a m and / a in L 2 (T") is straightforward from the theory 



of monotone operators; see Section 2.2 C 2,a regularity of / a>m follows from the elliptic regularity 
theory since / + adE m is a quasilinear uniformly elliptic operator; see [ fT5| . The Lipschitz 
bound on / a ,„ is obtained by a standard comparison argument. Indeed, since f ajn is a classical 



solution of ( |5.3b[ ) with right-hand side ef and the problem is translation-invariant, we have that 



f a ,m(- +z) -e HV/IL |z| is a subsolution of ( 53b ) with right-hand side s(f(- + z) - ||V/|L |z|) < e/ 



for any z € T". Therefore the comparison principle yields f a , m (x + z) - £ IIV/H^, |z| < f a , m (x) for 
all x,ze T", which implies the Lipschitz bound [[V/a^U^ < s IIV/H^. 

As we shall show, the approximation of the £ -resolvent problem ( |5.3a[ ) by the ^,,,-resolvent 



problems ( |5.3b| ) is a consequence of the Mosco convergence of £,„ to £. Here we present 
only the outline of the argument and we refer the reader to the book of Attouch pj for the 
definition and further details. Indeed, due to Proposition |5. [|5j Theorem 3.20] implies 
that E m Mosco-converges to E as m — » oo. Consequently, the Mosco convergence of E m to 
E implies the graph convergence of dE m to dE as m — > oo; see pi Theorem 3.66]. Finally, 
the graph convergence is equivalent to the resolvent convergence: we have, for fixed a > 0, 
f a ,m - » fa in L 2 as m — > oo; see [|5l Theorem 3.62]. We conclude by recalling the uniform 
Lipschitz bound ||V/a, m || < e HV/H^, and its consequence ||V/i ajm | < 2s HV/H^, /a, which yield 
the uniform convergence of f a<m and /z am as m — > 00. In particular, we recover the Lipschitz 
bound ||V/ fl |L<e||V/|L. 



By the classical theory of resolvent problems, f a — » e/ in L (T' ! ); see [21]. This yields, 
together with the uniform Lipschitz bound ( |5.4[ ), the uniform convergence. Additionally, conver- 
gence of h a is arestatement of Proposition |2.2| in light of Remark 2.10 



□ 



Proof of Stability Theorem 5.2 We show the statement for subsolutions. To simplify the nota- 
tion slightly, let us denote u = u. 

Clearly u m is an upper semi-continuous function. We must verify that (i) and (ii) of Defini- 
tion [U hold. 

(i) Let (p(x, t) = f(x) + g(t) be an admissible test function in a general position of radius 77 > at 
(x, i) with respect to u. By definition, / is the support function of some smooth pair + ) 
with Lipschitz constant L := WVfW^ < 0. Observe that the definition also ensures f(x - z) = 
for all |z| < 77. We set 6 := rj/4. 
Let us introduce the closed sets 

U = {u(; t) > u(x, f)}, Z = {/ < 0}, N p = U p n Z p for p > 0. 

The situation at time t is depicted in Figure |2j 

We will fix a number s e (0, 1) throughout the proof. The following observation is an 
important consequence of the general position of (p and u, but we will postpone its proof until 
after the end of the proof of stability. 

5.4 Lemma. Assume that <p is in a general position of radius n with respect to u. Then, by 
adding the term \t - f| to git) if necessary, we have 

u(x, t)-s inf f(x - z) - g(t) < u(x, t) - g(t) (5.6) 

|z|<<5 

whenever (x, t) e [N v x [i - 77, i + 77]) \ (N s x {t}), where 8 = rj/4. 
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z _ 
u 



Figure 2: Geometry of the graphs around the contact point (x, t), at time f, in the proof of 
Theorem 15 .21 

Now we recall the definitions of f a and f a m , together with h a and h a m , and their properties 
from Proposition [53] which we will use in the rest of the proof. 
For given a > and z, \z\ < 6, we define the set of maxima 

A z := arg max [u(x, t) - f a (x - z) - g(t)] . 

N 3 sX[t-26,t+2S] 



The uniform convergence f a =| ef and the strict ordering in Lemma 5.4 guarantee the existence 
of do > such that whenever a < oq then 

* A z c N 2S x [t -S,t + 6] for all | z | < 5. (5.7) 

We fix one such a < ao and find z, \z\ < 6, such that 

f a (x-z) = min/ a . (5.8) 

This choice will guarantee ( |5.1 1[ ) to hold. 

Again, due to the uniform convergence f ajn =4 f a as m — » oo, there is (x a , £ a ) e A z and a 
sequence (x fl)OT , ? ajOT ) (for a subsequence of m) of local maxima of 

(X, t) U m (x, t) - f a ^(x - z) - g(t). 

such that (x a>m , f a m ) — > (x a , £ a ) as m — > oo (along a subsequence). 
Since / fl)W e C 2 a and w m is a vise, subsolution of (|5.1[), we have 



,? (fa,m) fa,m(x a m Z), h a m ( y X a m z)) ^ 

— & ifa,m) fa,m(.X-a,m Z)j dE m (f a m )(x am Z))) ^ 0. 

Now we recall the Lipschitz bound ( |5.4| ), which yields that there exists p a such that, along a 
subsequence, V/ a>m (x Qim ) — > p fl , and |pj < eHV/H^, = eL. Consequently, we send m — > oo in 



(5.9) along a subsequence, and use the uniform convergence h am =| h a and the continuity of g' 



and F to conclude that 

*U) + F(p 8 A(jQ,-Z))<0. (5.10) 

To estimate h a (x a - z), we observe an important geometric consequence of the choice of z in 



(5.8 ), and the profile of u, which follows from the general position of ip with respect to u. Again, 



we postpone the proof until after the end of the proof of stability. 
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5.5 Lemma. We have 



h a (x a - z) < h a (x- z) = mmh a . (5.11) 

B 6 (x) 



Therefore, by elliptic ity of F and by ( |5.10| ), we arrive at 

g'(t a ) + F(p a ,minh a ) < g\t a ) + F(p a ,h a (x a - z)) < 0. (5.12) 

B 5 (x) 

In the next step, we send a — » along a subsequence. We choose the subsequence ai 
in such a way that p Ul — > /> with |p| < eL (this is possible by \p a \ < sL) and min^ h ai — > 
lim inf a _^ muig /z a as Z — > oo. Since also t ai — > f, we send Z — > oo and, with the help of 
continuity, we conclude from ( |5.12[ ) that 

g'(t) + F(p,liminf min/zj < 0. (5.13) 

It is time to use the knowledge of the limit h a — > -d°E(sf) in L 2 ; we remind the reader 
that e/ e D(dE), as was observed in Remark 2.10 Because B s (x) is contained in the facet, 
-d°E(sf) = A(Q_,Q + ) a.e. in B s (x) by Proposition |2~9j This yields 

lim inf minh a < ess inf A(Q_,Q + ). 

B 5 (x) Mx) 



One last time, the above inequality, together with the ellipticity of F and with ( |5.13| ), leads to 

g'{t) + F(p, ess inf A(Q_,^ + )) < 0. (5.14) 

As the last step, we recall that the reasoning is valid for any fixed s 6 (0, 1), yielding 
p = p(s). After sending s — » and realizing that \p(s)\ < sL — > 0, we recover the inequality 
( |5.14[ ) with p = 0. This finishes the proof for admissible faceted test functions. 

(ii) We can follow the standard stability argument in the theory of viscosity solutions. 

The proof for supersolutions is completely analogous. □ 



We still owe the reader a proof of Lemma 5.4 and 5.5 which both rely on the following 
geometrical lemma. 

5.6 Lemma. Suppose that <p is in general position of radius r] with respect to u at (x, t), where 
<p(x, t) = f{x) + g(t). Then 

u(x, t) < g(t) - g(t) + u(x, t) for (x, t) e TP X[t-rj,t + rj], 

and 

fix) > for x 6 U\ 

where Z := {f < 0} and U := {«(•, t) > u(x, i)}. 

Proof. By definition of general position, Definition 3.2[ 

u(x -z,t)- f(x) - git) < uix, t) - fix) - git) for all | z | < 77. 

We recall that fix) = 0, and, by definition, whenever y eZ 1 there exist x e Z and z, \z\ < rj, such 
that y = x - z. Similarly, if x e U n then there exists z, \z\ < r\, such that x - z & U. The result 
follows. □ 



28 



Proof of Lemma \5~4\ We shall use the notation from the proof of Theorem |5.2[ Let us recall 
that e e (0, 1). For any t e [t - 77, i + 77], we split the set N25 into two disjoint parts: 

• If x 6 N AS \ Z 6 c Z 4S = Zi, it is clear that inf^ fix - z) > and therefore 

u(x, t)-s inf f(x -z)- git) < u(x, t) - git) < u(x, t) - g(t), 

\z\<S 

where the second inequality follows from Lemma |5T6| and the fact that x e Z v . 

• On the other hand, if x e n Z s then inf| z |<5 f(x - z) < and the definition of general 
position yields 

u(x, t)-s inf f(x -z)- git) < u(x, t) - inf f(x - z) - git) < uix, t) - fix) - git). 

\z\<6 \z\<S 

One readily observes that 

(a) the first inequality is strict whenever inf^ fix - z) < (recall that s < 1); 

(b) the second inequality is strict whenever 

o t ^ f (by adding \t - f| to git) if necessary), 
o t = t , inf| z |<,5 fix - z) = and x i U. 



We recover Lemma 5.4 by collecting all these cases. 



□ 



Proof of Lemma 5.5 Recalling the definition of h a in ( |5.5| ), we want to show 

fa(Xa -Z)~ fix a - z) < f a (x - Z) - fix - z) = Hlin f a - f. 



The equality follows directly from ( |5.8| ), after realizing that fix, - z) = since x - z lies in the 
facet. 

To prove the inequality, we make the following two observations: 

(a) -f(x a - z) < = fix - z): indeed, since x a e N 2 s due to ( |5.7| ) and the choice of a, we 
observe that x a - z & N 3S c U v and therefore Lemma [53] applies. 

(b) f a ix a -z) < f a ix-z): as ix a , t a ) e A z , i.e ix a , t a ) is a point of maximum, and x e U C\Z c N is , 
we must have 

U(x a , t a ) - f a ix a ~Z)~ g(Q > ui%, t) - f a ix - z) - git)- 

A rearrangement of this inequality yields 

faiXa ~Z)< faix ~ z) + [u(x a , t a ) - g(t a ) - u(x, t) + git)] , 



and, because x a e Z v and thus Lemma 5.6 implies that the expression in the bracket is 
nonpositive. 



The inequalities in (a) and (b) together imply Lemma 5.5 



□ 
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5.2 Existence 



The standard consequence of the stability is the following existence theorem. 

5.7 Theorem (Existence). If F is degenerate elliptic and Uq e C(T"), there exists a unique 
solution u e C(T" x [0, oo)) of < |1 . ip with the initial data uq. Furthermore, ifuo e Lip(T n ) then 

IIVmC-, Olloo < IIVmoIIoo - 

The proof of the theorem will proceed in three steps: 1) due to the stability, by finding the 



solution u m of the problem ( |5.1 ) for all m > 1, we can find a subsolution u and a supersolution u 
of < | 1 . 1 1 ) ; 2) a barrier argument at t = shows that u and w have the correct initial data u ; and 
3) the comparison principle shows that u = u is the unique viscosity solution of and the 
Lipschitz estimate holds. 

Before giving a proof of the existence theorem, we construct barriers for step[2j 
Since the operator ( |5.2| ) degenerates at points where Vw = as m — > oo, it seems to be 
necessary to construct barriers that depend on m. We will use the Wulff functions for energy E m ; 



these were previously considered in the proof of stability for general equations of the type (5.1 ) 



in one-dimensional setting in [26 \. 
As we did before, let us set 

W(p) = \p\ , W m (p) = (\p\ 2 + m- 2 ) 1 ' 2 + m" 1 \p\ 2 . 

Given these convex functions, we define their convex conjugates (via the Legendre-Fenchel 
transformation) 

CW = sup [p ■ x - W m (p)] , x 6 R". (5.15) 

W*(x) is defined analogously. We recall that W(p) is the support function of 5i(0) and therefore 

f0 \x\ < 1, 

+oo Id > 1. 



Functions W* m are the Wulff shapes of the energies E m in the sense that the operator ( |5.2[ ) 
applied to W* m gives a constant for all x. We summarize the properties of the Wulff function in 
the following lemma. 

5.8 Lemma. The functions W* m are: (i) smooth functions on R"; (ii) strictly convex, radially 
symmetric; (Hi) increasing in m; (iv) W* n > -m~ l with equality only at x = 0; (v) — » 
uniformly on {\x\ < 1} as m — » oo; (vi) ~W* m — > +oo uniformly on {|x| > R} for any R > 1 as 
m — » oo; and finally, (vii) 

-d°E m (W* m ) = n. 

Proof. Since W m is smooth and D 2 W P > m~ l I, we have that W m is in fact uniformly convex 



and W m (p)/ \p\ -> oo as \p\ — > oo. Therefore the supremum in ( |5.15[ ) is attained at a unique 
point p such that x = V p W m (p), and p can be written as a smooth function of x using the 
implicit function theorem, (ii) is clear since W m strictly convex, radially symmetric; (iii) 



holds because W m are decreasing in m; (iv) follows from ( |5.15[ ) with p = 0; we observe that 
m~ l < W* n < in {|*| < 1}, which yields (v); (vi) can be derived from convexity of W^, ( |5.15[ ) 



and from the fact that W m (p) \ \p\ as m — > oo; finally, (vii) follows from the Fenchel identity 



x = VpW m (p) if and only if p = VW'(jc), which gives, using the formula in Proposition 5.1 ^b), 
-d°E m (W*J(x) = div [(y p W m )(VW* m (x))] =divx = n. □ 
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Since F depends on the gradient as well, while the gradient of W* m blows up as m — > 0, 
we consider rescaled cut-offs of these functions with linear growth at infinity. Given positive 
parameters A and q, where A will represent scaling and q the maximal slope, we introduce a 
smooth, nondecreasing cut-off function 



e q {s) 



\s \s\<q/2, 
I q sign s \s\ > 2q. 

Since W* n are smooth, radially symmetric, we define the rescaled modified Wulff functions on 



as 



w m (x; A, q)= f 6 q (Asx/ \x\)) ds. 

o 



5.9 Lemma. Let A and q be positive constants. Then w m (x;A, q) are smooth, convex, radially 
symmetric, nonnegative functions on R n , and they coincide with A~ l W* n (Ax) + (Am)~ l on the 
closed ball i\VW* n \(Ax) < q/l\. Moreover, \Vw(-,A,q)\ < q. 

Functions w,„ will be used to build barriers at t = for the existence argument. 

5.10 Lemma. Given positive constants A and q, there exist universal constants B = B(A, q, F, n) 
and mo = mo(A, q, F, n) such that, for all m > mo, 

<p m (x,t;A,B,q) := Bt + w m (x;A,q) 



is a (classical) supersolution of ( |5.1| ) in R" x [0, go), and -<p m (x,t;A,B,q) is a (classical) 
subsolution. 

Proof. We choose m large enough to guarantee that Vw m (x) = VW* n (Ax) in a neighborhood of 
x = 0, uniformly in m > m = 4/q. 



Let us now explain why m = 4/q is a reasonable choice. We first observe that Lemma 5.8 



implies that |VW^ (Ax) is an increasing function of |jc|, and that -(Am)~ l < A" 1 W* n (Ax) < 
in < A -1 }- From this, and the choice of m = 4/q, we can conclude that |VW*,| (Ax) < 
(Am)~ [ /(A~ l /2) = 2m" 1 < q/2 in l\x\ < A~ l /2}. Finally, Lemma 5.9 and the chain rule yield that 
Vw m (x) = VW* m (Ax) for all \x\ < A~ l /2 and m > m . 
To find B, we split R" into two sets: 

• if x such that |VW^(Ax)| < q/2: then Vw m (x) = VW* m (Ax) and thus we conclude that 

-d°E m (w m )(x) = An; 

• otherwise: then as observed above, we must have |jc| > A~ x /2 and |Vw m | > q/2, which 
yields \-dE m (w m )(x)\ < C for some constant C, C > An, independent of x and m > mo; 
this is a straightforward calculation, by evaluating the operator ( |5.2| ) applied to w m , which 
is radially symmetric, and observing that it can be bounded from above, independently of 
m and x, provided that |x| > A" 1 /2 and |VwJ > q/2. 

Finally, by construction, |Vw m | < q in R n . Since F is continuous, 

B := max \F(p,s)\ <oo. 

\p\<iM<c 

<p m are smooth and we have 

(<p m ) t + F(V0 m , -dE m (<t> m )) >B-B = 0. 



From this we see that (p m is a classical supersolution of ( |5.1[ ). Since -5£ m is odd, the same 
computation shows that -<p m is a subsolution. □ 
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Let us conclude by finishing the proof of the existence theorem. 



Proof of Theorem \5J\ Step 1. For every m > 1, there exists a unique solution u m of ( |5.1| ) with 
initial data uq. 

Since u Q is bounded, u m are uniformly bounded, by comparison principle, on sets T" x [0, T] 
for every T > 0. Therefore the following half-relaxed limits are bounded: 

u := ★-limsup u m , u := ★-liminf u m . 



Additionally, by Theorem 5.2 u is a supersolution and u is a subsolution of (jTTT|). Furthermore, 

u(x, 0) < Mq(x) < u(x, 0). 

Step 2. We need to show that in fact the above inequality is an equality. To achieve that, we 
will construct a class of radially symmetric sub- and supersolutions of ( |5.1| ) for all m, arbitrarily 
close to the initial data. 

We set K := 2 ||w |L < 00 • Recall that T" is compact, hence u has a uniform modulus of 
continuity. For any fixed 5 > 0, there exists r e (0, 1/8) such that 

|« (x) - Ko(£)l < <5 for all f,xe T, |x - f| < 4r. (5.16) 



We fix 5 > and set A = r _1 , g = 2Kr~ l . By Lemma 5.8 Vi), there exists m such that 
A" 1 W^(Ajc) > for all x e K", |x| > 2r, and m > m . We choose m > m so that the conclusion 
of Lemma 15 . 1 01 holds with a universal constant B. 
We claim that 

<f> m (x,0;A,B,q) = w m (x;A,q) > K for \x\ > 3r. (5.17) 

Indeed, by convexity, |VW^| (Ax) and \Vw a (-;A,q)\(x) are increasing functions of |x|. If 
|W;|(Ax) < q/2 at \x\ = 2r then w A (y;A,q) > w A (x;A,q) = A^W^Ax) + (Amy 1 > K 
for \x\ = 2r and \y\ > 2r. On the other hand, if |VW*|(Ax) > q/2 at |jc| = 2r, then 
w A (x;A,q) > |(|x| - 2r) for |x| > 2r, and the claim is proved since \i\x\ - 2r) > y = K 
for |jc| > 3r. 

For arbitrary x,£ 6 T", we set 

— ^ 

m (x, := -0 OT (x - £ A, 5, q) + i* (£) - ft 
$ (x, t) := m (jc - A, 5, #) + m (£) + ft 

— m 

where m is the "periodization" of (p m , i.e., 

<L(x, r; A, 5, o) = inf OT (x + z, t; A, 5, a), jeR",f>0. 

zeZ" 



We remark that 0,„ is still a viscosity supersolution of ( |5.1[ ), due to the stability properties of 
viscosity solutions. Since <p(x,0; A, > 0, 4r < 1/2, and ( |5.17[ ) and ( |5.16[ ) hold, we have 

0) < w < ^ (•, 0) in T" for £ e T", m > m . 

— m 

Therefore, by the parabolic comparison principle, 

~$ m < u m < in T x [0, oo), for £ e T", m > m . (5.18) 
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Recalling that 

f m (x, t) > -Bt + «o(£) -6- {Amy 1 for \x - £| < r, 

we have ★4iminf m _+ 00 m (£, 0) = « (£) - 8. 

We can similarly compute ★-limsup m „ >00 f (£, 0) = Uq(£) + 5 and we recover from < |5.18[ ) 

u Q (£) -S< «(£ 0) < w(£ 0) < h (£> + <5 for £ e T" . 
Since £ > was arbitrary, we have 0) = u{-, 0) = w - 



Stop 3. The comparison theorem 4.1 then yields u = u in T" x [0, oo), and it is the unique 
solution of (jTTTJ) with initial data u . 

To show the Lipschitz estimate, we use the standard comparison argument. Suppose that 
wo 6 Lip(T"), and set L = ||V«|L, and let u be the unique viscosity solution. Since the operator 
F does not depend on neither x nor u, we observe that u(x + h) + L \h\ is a viscosity supersolution 
for any h e T". Therefore, by comparison, 

u(x, t) < u{x + h,t) + L \h\ for x,he T", t > 0. 

The Lipschitz estimate follows. □ 



The proof of the existence theorem, Theorem 5.7 also gives the following approximation 
result. 



5.11 Theorem. Let uq e C(T") and let u m be the unique viscosity solutions of ( |5.1| ) on T"x(0, oo) 
for m > 1. Then u m converge locally uniformly on T' 1 x [0, oo) as m — > oo to unique viscosity 
solution u of < | 1 . 1 b vvzY/z initial data uq. 



Proof. It was shown in the proof of Theorem 5.7 that 



u = ★-limsup u m = ★-liminf u m , 
which is equivalent to the local uniform convergence. □ 



5.12 Remark. Theorem 5.11 implies that in the case of the standard total variation flow equation 



(1.4) the viscosity solutions defined in this paper coincide with the semigroup solutions given by 



the theory of monotone operators. 
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